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A model of weakly interacting hole quasiparticles is proposed to describe the nor- 
mal state of the high temperature superconductors. The effect of strong correlations 
q \ is contained in the dispersion relation of the holes, which is obtained accurately using 

a numerical technique and the t-J model on 12 x 12 clusters. Saddle-points gener- 
ated by many-body effects induce quasiparticle flat bands similar to those observed 
in recent angle-resolved photoemission experiments. The near degeneracy between 
momentum states along the X — Y direction is a hidden small energy scale in the 
problem, which produces a strongly temperature dependent Hall coefficient, as well 
as a large Fermi surface similar to those observed in the cuprates. 
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Many normal state properties of the high temperature superconductors are still puzzling. 
The in-plane d.c. resistivity, p a b, of the hole-doped cuprates is linear with temperature (T) 
when the hole concentration is optimal i.e. when the critical temperature (T c ) is the highest. 

The Hall coefficient, Rh, is strongly temperature dependent in the range between T c to 
T ~ 800K. At constant temperature, Rh changes sign as the hole density is increased 
away from the insulator parent compound. fHH Angle-resolved photoemission experiments 
(ARPES) suggest the presence of a large (electron-like) Fermi surface in some cuprates. @] 
This result is in contradiction with theoretical ideas at T=0 based on electronic models 
which tend to favor the presence of small hole-like Fermi surfaces at low hole density. || In 
addition, using ARPES techniques it has been recently reported || that an extended region 
of flat CuCVderived bands very near the Fermi energy exist for Bi2212, Bi2201, Y123 and 
Y124. This unexpected result cannot be explained within band structure calculations. || 

In this paper we discuss a model of hole carriers in an antiferromagnetic background 
that explains in a natural way many of the above mentioned anomalous properties of the 
cuprates. The main assumption is that the normal state regime of the cuprates can be 
modeled by a non-interacting gas of hole-quasiparticles. These carriers are strongly dressed 
by spin-waves. The influence of antiferromagnetism and strong correlations are contained in 
the special dispersion relation, e(k), which is obtained using an unbiased numerical method 
applied to a one band model description of the CuC>2 planes. Our approach is similar in spirit 
to that of Trugman where the dispersion relation of holes was calculated variationally 
using the tt'J model. One of the main results of the paper is that the hole dispersion relation 
contains flat bands at the X and Y points (i.e. k = (tt, 0), and (0, tt)), very similar to those 
described in the recent ARPES experiments. || The present approach belongs to the family 
of "narrow band" descriptions of the normal state of the cuprates, with flat bands caused by 
strong electronic correlations instead of band structure effects. In addition, e(k) contains a 
somewhat hidden energy scale, A, which is considerably smaller than the total quasiparticle 
bandwidth of the order of the antiferromagnetic exchange J. While this small scale was 
observed in previous calculations by Trugman, [0], Dagotto et al., f| and others, its physical 
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consequences have not been compared with recent experiments for the cuprates. We show 
that the strong temperature dependence of Rh and the Fermi surface (FS) induced by the 
presence of A produces results in good agreement with the data obtained by Hwang et al. 

1 

To calculate e(k), here we used the two dimensional t — J model. |J However, the in- 
tuitive ideas to be discussed below are not dependent on the specific model selected to 
represent the cuprates, as long as strong antiferromagnetic fluctuations are present in the 
ground state. Using the one band Hubbard model near half-filling, the results of our analysis 
are qualitatively the same. The t — J model Hamiltonian is defined as 

H = -t £[(1 - n w )4cj.(l - n w ) + h.c] + J^[Si-Sj - imm], (1) 

(ij> (ij> 

where the notation is standard. The energy of holes in this model has been extensively 
analyzed in the literature using a wide variety of numerical and analytical techniques. ||] 
Here e(k) is calculated using a Green Function Monte Carlo (GFMC) method ]TD[ which 



allows to obtain results on large clusters of 12 x 12 sites minimizing finite size effects. We 
have explored other techniques to carry out this calculation (Lanczos method on clusters 
from 16 to 26 sites, and the Born approximation at large spin S), and all of them produce 
results in good agreement with GFMC. 

In Fig.l, the numerically evaluated e(k) is shown for J/t = 0.4 along particular directions 
in the Brillouin zone. The total bandwidth, W, is severely reduced from that of a gas of 
non-interacting electrons due to the antiferromagnetic correlations. The minimum in the 
energy is obtained at the M point k = (ir/2, vr/2). These observations are in agreement with 
a wide variety of previous calculations, [|7|,[|] which have consistently obtained a bandwidth 
W ~ 2J and the energy minimum at M. However, the analysis below suggests that the 
strong anisotropy of e(k) is more important than the actual bandwidth for calculations 
carried out near room temperature. The close proximity in energy of momenta M and X 
introduces a small energy scale A = e(X) — e(M) in the problem in agreement with previous 
variational and numerical calculations. |7]|| Actually, all momenta belonging to the non- 



interacting 2D Fermi surface cosk x + cosk y = are very close in energy in the t — J model. 
In the range 0.3 < J/t < 0.7, A is approximately 15% to 20% of the total bandwidth. Since 
W ~ J, then important effects are expected at T ~ 300K i.e. room temperature. 

In order to calculate some observables, we have fit the numerical results of Fig.l using a 
general combination of trigonometric functions. The best fit corresponds to 

e(k) = -1.255 + 0.34 cosk x cosk y + 0.13(cos2k x + cos2k y ), (2) 

where higher order harmonics carry negligible weight. The main effective contribution 
to e(k) arises from hole hopping between sites belonging to the same sublattice, to avoid 
distorting the AF background. For completeness, we have also calculated numerically e(k) 
when a n.n.n. hopping term of strength t' along the diagonals of the square lattice is added 
to Eq.(l). Using t'/t = —0.15, [[llj the energy minimum is now located at X as shown 



in Fig.l, and the mass anisotropy between directions along and perpendicular to X-Y, is 



even higher than at t' = 0. |p~2|JT3[1 However, the actual position of the hole minimum is not 
crucial for our arguments below, but more important is the near energy degeneracy along 
the cosk x + cosk y = line. Thus, our conclusions are not affected by a nonzero t'. 

An important detail of Fig.l is the near flatness of the energy in the vicinity of k = (jr, 0). 
This feature is in striking agreement with ARPES results obtained by Dessau et al. || 
for Bi2212 (reproduced in Fig.2a). These authors remarked that such a flat region near 
(0, 7r), (it, 0) is observed in several high-Tc compounds and seems a universal property of 
the hole-doped cuprates. To compare the theoretically observed flat region of Fig.l and 
the experiments, we fixed the chemical potential of our model such that k = (tt, 0) would 
approximately be the Fermi momentum pp and to set the overall scale we used t = 0.4eV. 
Only momenta from T to the Fermi momentum are plotted since for momenta above pp, 
the data points of Fig.l represent poles in the hole spectral function that do not necessarily 
carry a large weight, and thus presumably they do not correspond to the large quasiparticle 
peaks observed in the ARPES experiment. This issue has been discussed at length in 
previous publications. 0] The qualitative agreement between the t — J model prediction 



and experiments is remarkable showing that the flatness of the k = (jr, 0) region in the 
cuprates may have a many-body origin. In the present calculation the presence of saddle- 
points near X and Y (also noticed in Ref. 0) seems responsible for the abnormal flatness 
of the hole dispersion. The position of the saddle points (SP) at J/t = 0.4 is approximately 
at k = (0.737T, 0) and their rotated ones. From Eq.(2), it can be shown that the location of 
the SP near the X point is a consequence of having the hole-band minimum at the M point. 
As expected, the density of states, DOS, presents a van-Hove singularity caused by the 
saddle-points (Fig. 2b) 0. The similarities between the present approach, and the van-Hove 
scenario widely discussed before in the literature are remarkable. |15| The main difference 
is that in the present calculation the generation of the saddle-point in the dispersion is 
through many — body effects. The consequences of this saddle-point in the d.c. resistivity 
and superconductivity will be discussed in a future publication. 

The presence of the small scale A produces interesting physical consequences. First, 
let us consider its influence on Fermi surface calculations. In Fig. 3a, the hole occupation 
number in momentum space, n^, is shown at 10% hole density, at several temperatures 
ranging from T = OK to T ~ 500K. It is clear that n^ at k « (tt/2, tt/2) rapidly decreases as 
the temperature increases. Reciprocally, the occupation near k = (tt, 0) and (0, ir) increases. 
Intuitively, this result is understandable since different regions in momentum space will 
start contributing appreciably to observables once the temperature is comparable to their 
energy with respect to the T=0 chemical potential. Thus, in our model the near degeneracy 
between M and X implies that a strong temperature dependence in many quantities should 
be expected. In Fig.3b,c the Fermi surface is shown as a function of temperature. At 
T=0 and a realistic hole density x = 0.10, the minimum of the hole band located at M 
implies the presence of hole-pockets in the Fermi surface (Fig. 3b). However, when the 
temperature becomes comparable to A all levels along the X-Y direction become equally 
populated and the hole-pockets are washed out. To characterize the Fermi surface at a fixed 
finite temperature, we decided to plot the lines in k-space where the Fermi distribution / is 
appreciably reduced from its maximum value, f max , which is always obtained at M. To fix 
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ideas, in Fig. 3c we plotted the line where the population of a given state has been reduced to 
only 10% of the maximum f max . At a temperature of the order or larger than A no vestige 
of the hole pockets remain, and the "Fermi-surface" is large and electron-like, as observed 
experimentally and discussed in the variational approach of Trugman. [|7] 

We have also investigated the effect of a small A on the one band Hubbard model. While 
at large U/t the results are quantitatively similar to the t-J model, at small or intermediate 
couplings the presence of A has more dramatic consequences. For example, at U/t = 4, 
A ~ 0.009t ~ O.OtMeV^ ~ 42K, i.e. approximately an order of magnitude smaller than for 



the t — J model at J/t = 0.4. |L3| Thus, the disappearance of the hole-pockets near half-filling 
will occur at an even lower temperature than for the t — J model. Although the Hubbard 
model at an intermediate coupling U/t = 4 may not be a good description of the high-Tc 
cuprates, our approach can explain interesting Monte Carlo numerical results reported in 
the literature. Working on a 16 x 16 cluster, T ~ 0.17t and electronic density (n) = 0.87, 
Moreo et al. have observed a large Fermi surface in a numerical simulation of the Hubbard 



model. |T(| However, our analysis shows that only at temperatures as low as T = 0.009t the 
near degeneracy between M and X can be resolved and thus for all practical purposes Monte 
Carlo simulations at currently accessible temperatures treat all the momenta along X-Y as 
degenerate in energy. Then, the absence of "hole-pockets" in the numerical analysis is now 



understandable, |L7j although other possible many holes effects that lead to a large Fermi 
surface can also contribute to the Monte Carlo simulation results. 

The value of n^ shown in Fig. 3a have important consequences for other experimentally 
measurable quantities. The change from a pocket-like to a large Fermi surface as the tem- 
perature is increased (at low hole density) suggests that Rh may also acquire a strong 
temperature dependence in the present model. To check this idea, Rh was calculated in the 
relaxation-time approximation as |T8| 



Rh = Vo £ k (vxfe-v y fe)(-j) 

N 2eh [E k vi(-§)P 1 J 

where N is the number of sites, Vo is the volume of the unit cell, and the rest of the notation 

6 



is standard. To avoid the complexities associated with interplane coupling between adja- 
cent Cu02 planes, here we compare results exclusively with La2_ x Sr x Cu04 at different Sr 
concentrations. In this compound, Vo = 3.8 x 3.8 x 6.6 x A 3 w 95 x (10 _8 cm) 3 . In Fig. 3d, 
the sign of the Hall coefficient is shown in the plane T-x. A sign-change occurs at low 
temperature as a function of x at a concentration x ~ 0.4, and near half-filling increasing 



the temperature at T ~ 400K. |L9[ Results for other values of J/t in the physically relevant 
regime are very similar. In Fig. 3d the region where the Fermi surface changes from hole to 
electron-like is shown. It is interesting to observe that there is an intermediate regime where 
Rh is positive while the FS is large, which is reminiscent of the experimentally observed be- 
havior in the cuprates. Fig.4a shows the hole density dependence of Rh at low temperature 
compared against recent experiments by Hwang et al. The agreement theory-experiment 
is excellent. Note that on the theoretical side, no other adjustable parameter is used besides 
J/t, and actually even the J/t dependence of the results is weak. Our simple approach 
correctly reproduces the order of magnitude and qualitative trends recently observed exper- 
imentally in LSCO. In Fig.4b, the temperature dependence of Rh is shown. At low density, 
the agreement theory-experiment is again very good. The shift of the experimental results 
may be caused in part by the use of polycrystals (single crystal measurements systematically 
present a smaller Rh II). More difficult to understand is the behavior at x = 0.22. While 
the overall scale of theory and experiment are similar, the qualitative trends are not. In 
our model, Rh changes sign not only as a function of hole density but also of temperature, 
while experimentally Rh decays slowly to zero as T is increased. Quite likely higher energy 
states not considered in our model are important to account for this experimental feature. 
Nevertheless, the overall agreement theory-experiment is remarkable. 

Summarizing, in this paper we presented a simple model of hole quasiparticles that 
accounts for several unusual properties of the normal state of the cuprates. The effect 
of strong correlations is contained in the dispersion relation of holes which is accurately 
evaluated using numerical techniques. Results are contrasted with previous calculations by 
Trugman where a variational method was used to calculate the hole dispersion relation. An 



overall nice agreement between the two approaches is found. A saddle-point appears in e(k), 
inducing a van-Hove singularity in the DOS of many-body origin. One of our main results 
is that this saddle-point produces a nearly flat hole band near the X point in agreement 
with recent ARPES experiments by Dessau et al. The near degeneracy between X and 
M momenta introduces a small energy scale in the problem, causing the Hall coefficient and 
Fermi surface to have a strong temperature dependence near room temperature compatible 
with recent experiments by Hwang et al. ||. 

We have recently received a paper by Bulut et al. [2(J where similar discussions for the 
quasiparticle dispersion are given. We thank these authors for sending us an advance copy of 
their paper. Conversations with D. J. Scalapino, A. Moreo, Z.-X. Shen, and R. Laughlin are 
also gratefully acknowledged. E. D. is supported by the Office of Naval Research under grant 
ONR N00014-93-0495. M. B. is supported by the Office of Naval Research under grant ONR 
N00014-J-92-1320. Numerical calculations were performed in part on a CM-5 supercomputer 
at NCSA, Univ. of Illinois, Urbana, Illinois (the technical contribution of Michael d'Mello 
(Thinking Machines Corporation) is acknowledged), and on an IBM workstation cluster at 
SCRI, Florida State University, Tallahassee, Florida. 
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FIGURES 

FIG. 1. Energy of a hole in the t — J model, e(k), vs momentum obtained with the GFMC 
method on a 12 x 12 lattice and J/t = 0.4 (in units of t). The numerical results are the full 
squares, and the solid line is to guide the eye. Note the flat region near the X point. The dashed 
results correspond to numerical results for the same lattice and J/t, but including a n.n.n. hopping 
term in the Hamiltonian, t' = — 0.15t. The error bars are not shown but typically they are ~ 0.02t 
at all momenta, with the exception of the the T and M points where they are ~ 0.20t. 

FIG. 2. (a) Energy vs momentum of the experimentally observed quasiparticle band obtained 
using ARPES techniques (from Ref. ||) (open circles). The dashed line is an extrapolation dis- 
cussed in ||. The theoretical predictions of the present paper are shown as full squares joined by a 
solid line. They correspond to J/t = 0.4, on a 12 x 12 cluster, (b) Density of states obtained from 
Eq.(2) (t — J model, J/t = 0.4) showing the van-Hove singularity. The unit of energy is t. 

FIG. 3. (a) Occupation number, n^, versus momenta at hole concentration x = 0.10, J/t = 0.4 
and several equidistant temperatures ranging from T = to T = O.lOt. The arrows indicate the 
direction in which T is growing. The rapid reduction in the hole population near M is apparent; 

(b) Fermi surface as defined in the text corresponding to J/t = 0.4, x = 0.1, and zero temperature; 

(c) Same as (b) but at T ~ 0.08t ~ 400K. The full squares denote the position of the saddle-points; 

(d) Sign of the Hall coefficient, Rh, in the plane temperature - hole density (T,x). The solid line 
separates hole-like (Rh > 0) from electron-like (Rh < 0) behavior, at J/t = 0.4. The dashed line 
separates the regime where the Fermi surface is small and hole-like from the region where the FS 
is large and electron-like. 



11 



FIG. 4. (a) Hall coefficient (Rh) (hi units of 10 _3 cm 3 /C) as a function of hole density (x). 
The solid line denotes results for the pure t — J model at J/t = 0.4 and T=0. The full squares 
are experimental results for La2- x Sr x Cu04 at T = 50K, taken from Hwang et al. (Ref. ||) ; (b) 
Rh against temperature at a fixed hole density. The upper solid curve is the theoretical result 
obtained at J/t = 0.4, and x ~ 0.12. The nearby full squares are experimental results from Hwang 
et al. (Ref. ||) corresponding to x = 0.125. The lower solid curve and accompanying full squares 
correspond to theoretical and experimental results, respectively, obtained at hole density x = 0.22. 
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